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Tracking of interplanetary spacecraft equipped with optical communications 
systems by using astrometric instruments is being investigated by JPL. Existing 
instruments are designed to work at night and, for bright sources, are limited by 
tropospheric errors. To provide full coverage of the solar system, astrometric track- 
ing instruments must either be capable of daytime operation or be space-based. The 
integration times necessary for the ground-based daytime photon statistical errors 
to reach a given accuracy level (5 to 50 nanoradians) have been computed for an 
ideal astrometric instrument. The required photon statistical integration times are 
found to be shorter than the tropospheric integration times for the ideal detector. 
Since the astrometric accuracy need not be limited by photon statistics even under 
daytime conditions, it may be fruitful to investigate instruments for daytime optical 
tracking. 


I. Introduction 

Several observables for spacecraft navigation based on 
a laser telemetry system are being investigated. One such 
observable is the difference in direction between the space- 
craft and a cataloged reference object in the same field 
of view of an astrometric telescope. The error for such 
a ground-based differential astrometric measurement in- 
cludes the directional error of the reference source, the 
photon statistical error, and the error induced by index of 
refraction fluctuations in the troposphere. The best condi- 
tions for making differential astrometric measurements are 
at night, when the low background light levels lead to small 
statistical errors, and at high angles of elevation above the 


horizon, where the troposphere errors are smallest due to 
the short path length through the atmosphere. 

Unfortunately, on any given day only a small portion 
of the solar system can be observed at high elevations dur- 
ing the night. In Section II, the daily available observa- 
tion times for the planets are examined for the cases in 
which observations are restricted to nighttime and limited 
in elevation. All the planets suffer seasons in which no 
nighttime viewing is possible, corresponding to periods of 
low Sun-Earth-Probe (SEP) angles. While nighttime an- 
gular position measurements might augment other observ- 
able types during a fraction of a mission, critical parts of 
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a mission may take place at low SEP angles due to launch 
criteria. The restricted nighttime visibility implies that 
ground-based astrometric observations of laser-equipped 
spacecraft may have to cope with the daytime sky back- 
ground. 

In Section III, the photon statistical error of an as- 
trometric measurement in the presence of a background is 
analyzed. In Section IV, some specific examples of the pho- 
ton statistical error are presented. In Section V, a simple 
model of the tropospheric error is discussed. The relative 
size of these errors is affected by the angular separation 
of the reference source and the spacecraft. A larger sepa- 
ration is more likely to allow the use of brighter reference 
objects and result in a smaller photon statistical error. A 
small separation would require the utilization of fainter ref- 
erence objects but would decrease the tropospheric error. 
The smallest usable field of view is set by the density of 
the catalog stars. It was expected that in the mid 1990s 
the Hipparchos mission would provide a catalog with an 
average of 2.5 stars per 1-degree by 1-degree field with ini- 
tial accuracy of 10 nrad per component and proper motion 
uncertainty of 10 nrad/year [1]. The average brightness of 
reference sources in a 1-degree by 1-degree field is magni- 
tude m v = 8 [2]. Using magnitude m v = 8 and 1 degree 
for the source-spacecraft separation, it is found that the 
integration times for daytime observations with 50 nrad 
accuracy are about 30 minutes. 

For the error models used here, the tropospheric error 
for such measurements is larger than the photon statistical 
error. However, the daytime photon signal-to-noise ratio 
is very unfavorable. Existing astrometric instruments, de- 
signed for nighttime operation, are not able to work in 
this regime. Limitations imposed by real detectors or by 
systematic effects such as sky brightness variation may be 
much larger than the photon statistical error. 


II, Nighttime Planetary Viewing 
Limitations 


Figures 1-7 are plots of the number of minutes per day 
that Venus, Mars, Jupiter, and Saturn are above 10 degrees 
or 30 degrees elevation with the sun below —15 degrees for 
the time span of 1990 to 2000. Mercury is almost never 
visible in the dark sky and is not plotted. Venus is al- 
ways at low elevation, and nighttime astrometric tracking 
is possible less than half of each year. The outer planets 
are unavailable for 25 percent or more of each year. These 
figures suggest that visibility limitations will severely af- 
fect the utility of astrometric tracking if it is limited to the 
nighttime hours. 


111. Photon Statistics of an 
Astrometric Telescope 

This section presents an estimate of the accuracy lim- 
itation placed by photon statistics on an astrometric tele- 
scope in the daytime. This category includes Ronchi-ruled 
telescopes and charge-coupled device (CCD) instruments. 
An ideal detector is able to record the position on the fo- 
cal plane for each detected photon. Any real instrument 
will have larger photon statistical errors than this ideal in- 
strument. CCD instruments are capable of dividing the 
image of a point source into many pixels and can approach 
the ideal detector scheme. However, limits on the size of 
CCD arrays limit the fields of view of such instruments. 
In a Ronchi telescope, a moving ruling is used to modulate 
the light incident on the detector with position informa- 
tion derived from the detected modulation [3]. For present 
designs, the field of view must be larger than the image of 
the star. This makes Ronchi telescopes more susceptible 
to background light problems than CCD instruments. 

The ideal astrometric telescope tracks a source for a 
time T and records the plane-of-sky coordinates (£, ff) = ^ 
for each photon detected. The troposphere and instrument 
resolution cause the photons to be smeared in the plane of 
the sky. For simplicity, the photon spatial distribution 
about the true source direction is assumed to have the 
Gaussian form 


To demonstrate what fraction of a mission trajectory ^ ^ ^ ST (fl — fi 5 ) 2 

would be inaccessible to nighttime astrometric measure- = W CXP 2a* (1) 

ments, the number of minutes per day that each planet is L 

visible from Goldstone in the night sky has been computed. 

Here “night” has been defined as the time that the sun is where B is the number ofjbackground photons per unit 

below —15 degrees elevation. This value for sun elevation solid angle per unit time, Q s is the true source direction, 

was chosen to correspond roughly to the time of astro- and S is the number of signal photons per unit time. This 

nomical darkness. Since the astrometric error depends on distribution is a reasonable approximation to a smeared 

elevation angle through the elevation dependence of the point source in a uniform background and has been suc- 

troposphere, two different minimum elevation cutoffs for cessfully used in fitting high-precision small-field images 

the planets were included. [4]. The width of the Gaussian is determined by the ap- 
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parent diameter of a point source. This apparent diame- 
ter as determined by the turbulence of the atmosphere is 
commonly called the seeing angle. The full width at half 
maximum of the Gaussian (2.35 <j) is here set equal to the 
seeing angle. 

Measurement of the source direction consists of 
recording the arrival of n photons and the plane-of-sky 
coordinates for each detected photon. The maximum- 
likelihood method can be used to estimate the source di- 
rection. The likelihood function L(fi a ) is the probability 
of recording n photons with the set of arrival directions 
{Qk} given an assumed source direction 

The photons are distributed in time according to Pois- 
son statistics. The expected number of detected photons 
is given by 

N = J /(Q)d 2 n (2) 


where < 7 2 is the variance of the source direction estimate 
in the £ direction, and the expectation value is the average 
taken over all possible numbers of detected photons and 
sets of directions. Since the assumed photon distribution 
is symmetric, the error is the same for the rj direction. 
With some work, the error expression of Eq. (6) can be 
applied to the likelihood function given in Eq. (5) to give 
the error in the source direction measurement as 



where the function f(a ) is given by 


/(«) = 



a + e 


dx 


( 8 ) 


In the limit of no background ( B = 0), a = 0, and /( 0) — 
1; the directional error is then given by 


where the integral is taken over the field of view of the 
telescope. The probability of detecting n photons is then 
given, in [5] as 


P{n) 


n! 


(3) 


The (normalized) probability of any one photon arriving 
with direction flfc given that the source is at direction Cl a 
is 


p{ 6*100 


nth - 0 ,) 

N 


(4) 


Assuming that the distribution for each photon-arrival di- 
rection is independent, the probability of finding n photons 
with the set of photon directions {Qfc } is 



J(Qib - fi„) 

N 


(5) 


The maximum-likelihood estimate of the source di- 
rection is the assumed source direction that maximizes the 
likelihood function. In the limit N > 1 the error asso- 
ciated with this estimate is given by the Konig-Kramer 
bound [6] 


d 2 log L(fi 0 ) 

da 


( 6 ) 


= 


a 

VST 


(9) 


Rewriting Eq. (9) gives the integration time to reach a 
specified accuracy as 


T = 



( 10 ) 


In the limit of high background, the position error is given 

by 


<7 /8xcr 2 R 

“ \/ST V S~ 


( 11 ) 


Rewriting Eq. (11) gives 



The effective signal-to-noise ratio for determining the 
source direction is 5/(87 T<t 2 5) where the combination 
87 t<j 2 B is the effective background rate. 


IV. Examples of Photon Statistical 
Integration Times 

In this section, the integration times needed to reach 
a given accuracy for a reference star and for a spacecraft 
at 10 AU with a 2-W laser are given. The parameters as- 
sumed for the detector, the spacecraft, and the background 
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are summarized in Table 1. The derived photon rates and The background photon rate is given by 

integration times are listed in Table 2. 

The photon rate from a star of visual magnitude m v (8na 2 )B = (8 n<r 2 ) I ( ) AA j d 2 r) ro r) ra r)jT] d (15) 

is given by [2, 7] as approximately \' lc / 4 


5=10 ° 4m *' 7 45 --- --- (A) AA j d^aVrartroV/Vd 

jim m l \hc } 4 

(13) 

where A is the central wavelength, h is Planck’s constant, 
c is the speed of light, AA is the wavelength pass band, 
d r is the telescope diameter, r\ a is the atmosphere trans- 
mission factor, 7] ra is the receiver obscuration factor, T] ro is 
the receiver optics efficiency, 77 / is the narrow-band filter 
efficiency, and r)d is the detector quantum efficiency. Since 
refraction will cause image smearing if the wave band is 
too wide, a value AA = 0.1 // m, centered at A = 0.532 fi m, 
will be used below. For the receiver efficiencies given in 
Table 1 and for a star of magnitude m v — 8 the signal rate 
is 6.8 x 10 5 photons/second. 

The photon rate received from the spacecraft is given 
by [7] as 


where I is the background spectral irradiance. The 
daytime background spectral irradiance depends on the 
weather, the sun elevation, and the Sun-Earth-Probe 
angle among other factors. The value for I of 
100 W/(^m m 2 steradian) at 0.532 /im [12] will be used 
in the following example. Kerr gives values ranging from 
30 W/(/im m 2 steradian) at 90 degrees from the sun to 
500 W//un m 2 steradian at 10 degrees to the sun [8]. Us- 
ing the receiver values of Table 1 and 2-arcsec seeing, the 
daytime background rate is 2.6 x 10 9 photons/second. Us- 
ing a narrower filter will not improve the signal- to- noise 
ratio for the reference star since the signal photon rate 
depends on the bandwidth in the same way as the back- 
ground; narrowing the filter bandwidth reduces the signal 
rate and increases the needed integration time. However, 
the laser signal is narrow band and narrowing the filter 
bandwidth improves the spacecraft signal- to- noise ratio. 


S = Pt 




X VtaVtoVtpVa^lra'nro 


(14) 


where Pj is the transmitted laser power, A is the wave- 
length, d t is the transmitter objective diameter, r is the 
spacecraft-receiver distance, rj ta is the transmitter obscu- 
ration factor, r]to is the transmitter optics efficiency, and 
Tj tp is the transmitter pointing efficiency. The parameters 
used here are taken from an example by Kerr [8] and listed 
in Table 1. These factors combined with the nominal re- 
ceiver used above yield a detected photon rate S of 2.6 x 10 4 
photons/second for the spacecraft. 

A critical parameter for the tracking telescope is the 
atmospheric seeing. At many sites, the daytime seeing is 
worse than the nighttime seeing by a factor of 5 to 10 [9, 
10]. However, at some solar observatory sites, the seeing 
is 5 to 10 /xrad (~1 to 2 arcseconds) both day and night. 
At the Sacramento Peak site, the daytime seeing angle is 
reported as being less than 10 /irad 80 percent of the time 
and better than 5 /irad 50 percent of the time during the 
day [11]. Since 10 /irad is similar to other reported daytime 
seeing values (see [11] and references therein), 10 /irad full 
width at half maximum will be taken as a nominal value. 
This corresponds to a value a — 4.3 j/rad for the photon 
distribution. 


Table 2 presents the integration times needed to reach 
50-nrad or 5-nrad accuracy for the reference star and the 
spacecraft. Several different filter bandwidth options have 
been used. For the case of the 0.03-nm filter, the filter 
transmission efficiency is reduced to 0.4 [8]. The inte- 
gration times for the photon statistics error to reach the 
5-nrad level are not prohibitively long provided that the 
spacecraft and star are separately filtered. However, sev- 
eral factors may combine to lengthen the integration times. 
The photon rates for the reference star plus background are 
very high for photon counting devices. Detector dead time 
and noise will be significant effects. A Ronchi telescope in- 
tegrates the background over a fixed field of view. Unless 
the field of view is limited to the size of the source image 
and precisely positioned, the photon statistics error will 
be much worse for a Ronchi telescope. In any case, the 
ruling modulation reduces the photon flux by a factor of 
2, correspondingly increasing the integration time. If each 
coordinate is measured separately, the times will double 
again. 

Effects other than counting problems are important. 
Since the integration time increases as the fourth power 
of the seeing, the site selection is critical. For the ref- 
erence star, the signal-to-noise ratio may be worse than 
10 -3 . This may pose unrealistic requirements on the dy- 
namic range and linearity of the detector. There is also 
the possibility that variations in the background intensity 
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over the period of integration could degrade the accuracy 
of the star position. 

V. Tropospheric Integration Times 

The photon statistics contribute an angular position 
error for the spacecraft and the reference object. The error 
in the angular difference between the spacecraft and refer- 
ence object is also affected by variations in angle of arrival 
imposed by the troposphere. The integration time require- 
ments imposed by the troposphere may be computed from 
a result by Lindegren [13]. Lindegren’s result is derived 
from a frozen turbulence model of the atmosphere with 
a power law structure function. The turbulence causes 
angle of arrival variation as the frozen turbulence moves 
with the wind velocity. The variations in angle of arrival 
from two point sources are computed, given some assump- 
tions about the atmosphere structure and averaging over 
wind direction. Lindegren’s paper compares his model to 
several experimental results, including stellar position and 
solar diameter measurements, with reasonable agreement. 

The expression for the error <7# in the difference angle 
0 between two sources is 

a e = 6.3 x Kr 6 0*T~* (16) 

where <r$ and 9 are given in radians and the integration 
time T in seconds. Solving Eq. (16) for the integration 
time necessary for a given angular accuracy <t $ gives 

T= 4.0 x lO” n 0*(T- 2 (17) 


This expression applies near zenith and with about 1- 
arcsec seeing. Results are expected to be worse for lower 
elevation angles and worse seeing. 

For two objects 1 degree apart, this implies an in- 
tegration time of 59 hours to reach 5-nrad accuracy and 
35 minutes to reach 50-nrad accuracy. These times are 
long compared to those imposed by the photon statistics 
for a magnitude 8 star. There is a possible trade-off by 
narrowing the field of view and utilizing fainter reference 
stars. However, this implies a larger catalog effort. By 
including more reference objects in the field, there is the 
possibility of reducing the troposphere error. 


VI. Discussion 

The photon statistics do not rule out the operation 
of astrometric telescopes for daytime optical tracking of 
spacecraft. Given the simple assumptions used in Sec- 
tion III, the tropospheric error dominates the photon sta- 
tistical error for a magnitude m v = 8 star and 1 -degree 
source-spacecraft separation. However, existing instru- 
ments are not capable of operating in the daytime because 
of the extremely poor signal- to- noise ratio and large pho- 
ton fluxes. Finding a suitable detection scheme for day- 
time operation will be a challenge. Since the integration 
times required are many minutes, systematic effects and 
variations in the background level may wash out the posi- 
tion signal. It would be useful to try differential position 
measurement in the daytime at a solar observatory with 
a CCD camera to begin an investigation into systematic 
background effects. 
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Fig. 1. Plot of the number of minutes per day (24 hours) for 
which Venus is visible above 10 degrees in the nighttime sky. 



YEAR 

Fig. 2. Plot of the number of minutes per day (24 hours) for 
which Mars is visible above 10 degrees In the nighttime sky. 



Fig. 3. Plot of the number of minutes per day (24 hours) for 
which Mars is visible above 30 degrees in the nighttime sky. 
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Fig. 4. Plot of the number of minutes per day (24 hours) for 
which Jupiter is visible above 10 degrees In the nighttime sky. 


96 







NUMBER OF MINUTES PER DAY NUMBER OF MINUTES PER DAY 



Fig. 6. Plot of the number of minutes per day (24 hours) for 
which Saturn is visible above 10 degrees in the nighttime sky. 
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